Pharmacology Research & Perspectives

Background: Diabetes mellitus is an increasing public health challenge in Sierra Leone, where access to diagnosis, treatment, and long-term care remains limited. Traditional medicine continues to play a significant role in disease management; however, ethnobotanical knowledge related to diabetes remains insufficiently documented. Methods: A cross-sectional ethnobotanical survey was conducted among 40 informants, including traditional healers, herbalists, and knowledgeable community members in Waterloo, Pendembu, and Bo. Data were collected using structured questionnaires administered via Kobo Toolbox and paper-based tools. Information on medicinal plants, plant parts used, preparation methods, routes of administration, and knowledge transmission pathways was obtained. Quantitative ethnobotanical indices, including Frequency of Citation (FC), Relative Frequency of Citation (RFC), and Informant Consensus Factor (ICF), were calculated. Results: A total of 21 medicinal plant species were documented. The most frequently cited species were Moringa oleifera (FC = 9; RFC = 0.225), Vernonia amygdalina (FC = 7; RFC = 0.175), and both Cassia siberiana and Telfairia occidentalis (FC = 6; RFC = 0.150). Leaves were the most commonly utilized plant part (40.9%), and decoction was the predominant preparation method (76.2%), with oral administration accounting for 95.2% of use. The Informant Consensus Factor (ICF = 0.69) indicated a relatively high level of agreement among informants. Knowledge was primarily transmitted through apprenticeship and inherited family practices. Conclusion: Traditional medicinal plants remain an important component of diabetes management in Sierra Leone. The high level of consensus among informants and the repeated citation of specific plant species suggest structured and culturally validated therapeutic practices. The findings provide a foundation for future phytochemical and pharmacological investigations and highlight the need for documentation, preservation, and sustainable utilization of ethnobotanical knowledge.

Background. Compounded tirzepatide is prescribed at scale as a cheaper substitute for branded Mounjaro and Zepbound, yet the cost case is almost always built by setting one compounded price against one branded list price. That framing ignores the question that actually decides the answer: cheaper than which branded price the patient can reach. Branded tirzepatide is now sold at sharply different tiers, namely insurance copay (often $25-$150/month), LillyDirect Self Pay ($299-$449/month), and retail cash price ($1,000-$1,200/month). Whether compounded saves money turns entirely on which of these a given patient faces. A second open question is whether the two formulations even produce comparable effectiveness, since observed differences may reflect selection on insurance, baseline characteristics, and adherence rather than the drug. Methods. We conducted a retrospective cohort study of tirzepatide users in the Mochi Health telehealth program, classified by formulation from their refills as branded-only (Mounjaro/Zepbound; 6,238), compounded-only (71,683), or switchers (4,996); switchers were excluded from the formulation contrast. Among single-formulation patients with a documented six-month weight observation, the analytic cohort was 7,271 (869 branded, 6,402 compounded). The primary outcome was six-month percent body weight loss; the secondary outcome was >=10% response. We used 1:1 nearest-neighbor propensity-score matching (0.25 SD caliper) on baseline covariates only - age, sex, baseline BMI, baseline weight, comorbid diabetes, hypertension, dyslipidemia, prior bariatric surgery, and self-reported insurance coverage - deliberately excluding post-treatment variables such as adherence and time in program, which are mediators of the formulation effect. We pre-specified an equivalence margin of +/-2 percentage points on mean loss and tested equivalence with two one-sided tests (TOST). A directed acyclic graph (DAG) makes the identifying assumptions explicit; metformin use could not be reliably ascertained and is treated as an unmeasured confounder. The cost comparison reports the savings or premium of compounded versus branded under five branded price scenarios: retail list, LillyDirect Self Pay (two dose tiers), and insurance copay (typical and low end). It is a cost comparison (cost-minimization under demonstrated similar effectiveness), not a formal cost-effectiveness analysis: we computed no ICER, QALY, or discounting. Results. Branded and compounded patients had similar outcomes even before adjustment (mean loss 11.7% vs 11.5%; >=10% response 60.9% vs 58.8%). The largest baseline difference between the groups was insurance coverage (branded patients far more likely insured; standardized mean difference 0.67), which matching balanced to 0.01. After 1:1 matching (718 pairs, all |SMD| < 0.04), mean loss was 11.4% vs 11.4% (difference +0.08 pp, 95% CI -0.70 to +0.80) and >=10% response 59.3% vs 57.2% (difference +2.1 pp, 95% CI -3.1 to +7.1). The two formulations were statistically equivalent within the pre-specified +/-2 pp margin (TOST p < 0.001). Cost depends on the branded scenario: compounded saves $6,000 over six months versus retail list price, $1,494 versus LillyDirect maintenance-dose (5-15 mg) Self Pay, and $594 over a low-dose (2.5 mg) LillyDirect prescription, while it costs $300 more than branded under a typical insurance copay ($150/month) and is more expensive still at lower copays (savings turn negative below $200/month). Conclusions. Branded and compounded tirzepatide were statistically equivalent in six-month effectiveness within a pre-specified +/-2 pp margin, so the choice between them is essentially a cost decision - and that cost advantage is real but conditional on the branded price the patient can access. It is large against retail list price and shrinks to zero or reverses against LillyDirect Self Pay or a low insurance copay. Whether compounded is the lower-cost choice for an individual patient is, therefore, a question about which price tier that patient faces.

Background Liver transplantation is a complex surgical procedure featuring prolonged operative time and extensive surgical trauma, which results in a high anesthetic risk. Especially during anesthesia induction and the anhepatic-to-reperfusion phases, where marked hemodynamic fluctuations may readily lead to malignant cardiovascular events. Liver transplantation is associated with numerous postoperative complications, including pulmonary complications, postoperative delirium, acute kidney injury, and delayed emergence, all of which may adversely affect patient prognosis. Studies on remimazolam besylate (hereafter referred to as remimazolam) have suggested that it has minimal impact on patient hemodynamics. In theory, this renders remimazolam an ideal sedative agent for liver transplantation. This study aims to verify the hypothesis that the use of remimazolam during liver transplantation can reduce the incidence of postreperfusion syndrome (PRS). In addition, we further explored the effects of remimazolam on postoperative complications in liver transplant recipients, particularly focusing on perioperative liver and kidney function, pulmonary complications, and postoperative delirium. Methods This study is a prospective randomized controlled trial. 120 participants aged 18-60 years who are scheduled to undergo liver transplantation under general anesthesia will be enrolled. In the intervention group, remimazolam besylate will be used for anesthesia induction and maintenance at doses of 0.2-0.4 mg/kg and 1-3 mg/kg/h until the end of surgery. In the control group, propofol will be used for anesthesia induction and maintenance at doses of 1-2 mg/kg and 4-12 mg/kg/h until the end of surgery. In both groups, anesthetic drug doses or sevoflurane administration for intravenous-inhalational combined anesthesia will be adjusted based on vital signs and BIS values. All other anesthetic medications will be conducted according to the anesthetist's preference and remain consistent. Discussion This trial will investigate whether remimazolam besylate can be used during liver transplantation to reduce the incidence of postreperfusion syndrome. It will also examine whether the drug provides potential benefits regarding perioperative complications such as postoperative delirium and acute kidney injury. Trial Registration: Chinese Clinical Trial Registry,ChiCTR2500095774, registered on January 13, 2025 Keywords: liver transplantation, remimazolam, postreperfusion syndrome, complications

Background. GLP-1 receptor agonist trials are tightly controlled: standardized titration, intensive dietary counseling, frequent in-person follow-up, and rigorous exclusion criteria. The real world is none of those things. In a U.S. telehealth GLP-1 program, diet engagement, exercise, medication choice, dose timing, and out-of-pocket cost vary substantially from patient to patient. Whether trial-level efficacy translates into the outcomes a patient and clinician will actually see is an open question, and the answer matters, because telehealth is now where most GLP-1 prescribing happens. Methods. We conducted a retrospective cohort study of 13,507 adults who used a single GLP-1 agent (tirzepatide or semaglutide) through the Mochi Health telehealth obesity program and had a documented six-month weight observation. The primary outcome was achievement of >=10% total body weight loss at six months. To address selection bias in the tirzepatide-semaglutide comparison, we used 1:1 nearest-neighbor propensity-score matching on age, sex, baseline BMI, baseline weight, and comorbid diabetes, hypertension, dyslipidemia, and prior bariatric surgery (recorded at intake), with a 0.25 SD caliper on the propensity logit. We drew a directed acyclic graph (DAG) with a clinical co-author to make the identifying assumptions explicit and to mark where unobserved variables (insurance, socioeconomic status, concomitant medications such as metformin) limit causal interpretation. We report multivariable predictors via logistic regression, compute an E-value for the matched contrast, and benchmark our point estimates against landmark RCT outcomes. Results. Overall, 59.1% of patients achieved >=10% loss at six months, with mean loss of 11.5% (median 11.3%). Threshold attainment was 86.6% at >=5%, 59.1% at >=10%, 27.5% at >=15%, and 9.1% at >=20%. The unadjusted tirzepatide-semaglutide response gap was +16.0 percentage points (68.8% vs 52.8%); after 1:1 propensity-score matching (3,480 pairs, all post-match |SMD| < 0.05) the gap was +18.1 percentage points (69.6% vs 51.6%, 95% CI +15.9 to +20.3). Matching on the measured covariates did not attenuate the advantage, indicating that selection on those characteristics does not explain it; the matched risk ratio was 1.35 (E-value 2.04). The gap was unchanged when a self-reported insurance indicator was added to the matching (+18.4 pp) and remained large (+14.2 pp) within patients who reached a therapeutic dose. Multivariable predictors of response were tirzepatide (OR 2.10, 1.95-2.26), female sex (OR 1.37, 1.20-1.56), and prior bariatric surgery (OR 1.36, 1.18-1.57); response was lower with comorbid diabetes (OR 0.84, 0.77-0.92) and, modestly, with higher baseline BMI per unit (OR 0.98, 0.97-0.99). Response varied by baseline BMI, from 58.0% in overweight patients (BMI <30) and a peak of 63.5% in Obese I to 52.4% in Obese III. Conclusions. Real-world response to GLP-1 therapy in a telehealth setting is meaningfully attenuated from RCT benchmarks but remains clinically substantial: roughly three in five patients reach the 10% threshold. The tirzepatide advantage over semaglutide is large and, notably, does not shrink under propensity-score matching on measured confounders, so it is not an artifact of the observed selection variables; an unmeasured confounder would need a risk-ratio association of about 2.0 with both drug choice and response to explain it away (E-value 2.04). The findings are observational, conditional on the DAG's identifying assumptions, and unmeasured confounders (insurance, socioeconomic status, concomitant medications) remain possible.

Polysubstance use is rising and linked to heightened overdose rates and increased treatment challenges, further exacerbated by increasing detection of adulterants (e.g., xylazine) in the street drug supply. Harm reduction groups provide sterile syringes in exchange for used ones, creating a unique opportunity to characterize prevalent polysubstance combinations and inform translational and preclinical research We analyzed residues from used syringes (N=3,168) obtained from several harm reduction organizations in Jefferson County, KY (Jan-Dec 2025) for the presence of substances using gas chromatography mass spectrometry (GC-MS). We classified compounds as adulterants (e.g., diphenhydramine [DPH]/Benadryl), byproducts/precursors of synthesis (e.g., 4-ANPP), and recreational drugs (e.g., meth). We excluded byproducts/precursors and determined the most frequent substance and pairs/trios containing one or more recreational substance. Results. Of 3,168 syringes, 2,522 (79.61%) tested positive for substances. Out of those positive, the top recreational substances were meth (n=1,387; 54.99%), fentanyl (n=1,220; 48.37%), and heroin (n=653; 25.89%). Top adulterants were DPH (n=1021; 40.48%), dimethyl sulfone (n=749; 29.69%), and lidocaine (n=736; 29.18%). The most common pairs were DPH+fentanyl (n=670; 26.57%), lidocaine+fentanyl (n=659; 26.13%), dimethyl sulfone+meth (n=621; 24.62%), and fentanyl+heroin (n=484; 19.19%). The most common trios were DPH+lidocaine+fentanyl (n=369; 14.63%), DPH+fentanyl+heroin (n=327; 12.97%), lidocaine+fentanyl+heroin (n=297; 11.77%), diphenhydramine+xylazine+fentanyl (n=273; 10.82%), and meth+lidocaine+fentanyl (n=262; 10.39%). Our findings highlight evolving patterns of multiple-opioid and opioid-stimulant polysubstance use, generating insights that can be rapidly applied to strengthen clinical, preclinical, and translational polysubstance research. These insights allow for investigations into biobehavioral mechanisms and consequences of emerging use patterns, accelerating development of novel therapeutics.

Background: We sought to determine the changing prevalence, incidence, and temporal trends in real-world, recorded diagnoses of metabolic dysfunction-associated steatotic liver disease (MASLD) and assess availability of fibrosis risk stratification following awareness campaigns and guideline updates over the last decade. Methods: This population-based cohort study identified MASLD diagnoses made between 2003 to 2022 in the UK primary-care Clinical Practice Research Datalink (CPRD) to estimate prevalence and incidence. A nested case-control analysis, utilising 1:4 age-, sex-, and general practice-matched controls, assessed clinical characteristics, availability of Fibrosis-4 (Fib-4) components, and its temporal trend pre- and post-2015. Findings: 11.7 million individuals were active in CPRD in 2022. 365,797 comprised the study cohort of people with a MASLD diagnosis (matched to 1,460,288 controls). From 2012-2022, recorded MASLD prevalence rose from 0.52% [N=51,028] to 2.42% [N=283,762] (p<0.001); recorded incidence doubled from 1.60 to 3.31 per 1000 person-years (p<0.001). People with MASLD diagnosis had a higher prevalence of type 2 diabetes (21.0% [N=76,640] vs 7.7% [N=112,812]) and hypertension (35.3% [N=129,156] vs 18.7% [N=273,502]). People of South Asian ethnicity were overrepresented in MASLD cohort but had the lowest availability of Fib-4 components (14.6% [N=4,467]; adjusted odds ratio 0.67, 95% CI: 0.65-0.70, vs White). Overall, Fib-4 availability increased pre- to post-2015 (4.3% [N=4,945] to 22.8% [N=56,634]). Among those with a calculable score, fewer South Asian individuals had indeterminate/high risk (18.6% [N=833] vs 35.3% [N=15,115] in White individuals, p<0.001). Interpretation: Recorded MASLD prevalence has increased 5-fold in a decade, yet a diagnostic gap persists. Fibrosis risk stratification has improved, but remains low and is potentially inequitable for people of South Asian ethnicity. Funding: Merck Sharp & Dohme LLC, a subsidiary of Merck & Co., Inc., Rahway, NJ, USA; Barts Charity. Keywords: Epidemiology, Real-world data, Real-world evidence, MASLD, Primary Care

Background. When a GLP-1 patient stops responding, should the clinician push the dose or change the drug? Observational answers conflate two distinct sources of confounding. Most early-week "escalations" in real-world data are FDA-mandated titration steps rather than deliberate clinical decisions, and patients who deviate do so for reasons we cannot observe. Semaglutide and tirzepatide are also not equivalent milligram-for-milligram, so naive class-switch comparisons mix mechanism and dose. We resolve both by restricting to post-titration patients and reclassifying treatments under the Whitley 2023 dose-equivalence framework. Methods. From 68,969 telehealth GLP-1 patients we built a post-titration cohort. Each patient's index time is the day they completed at least four weeks at therapeutic dose (Whitley tier 3 or higher: semaglutide 1.0 mg or tirzepatide 5 mg). Confirmed slow response is less than 5% total weight loss at the index, consistent with FDA weight-management drug-development guidance and AACE/ACE criteria. We compared four post-index strategies against continuing the current regimen: within-class dose escalation, equipotent class switch (a Whitley tier change of 1 or fewer), and class switch with potency increase. Direction-specific analyses split switches into semaglutide-to-tirzepatide and tirzepatide-to-semaglutide arms. Outcomes were percent weight loss at 12 and 24 weeks post-index. We estimated effects six ways: propensity-score matching; IPTW with linear and gradient-boosted propensities; the g-formula with linear and gradient-boosted outcome models; and AIPW, the doubly-robust estimator we use as the tiebreaker. Two-layer inverse probability of censoring weighting addressed strategy adherence and outcome ascertainment. We computed E-values, ran a negative-control specification, and stratified by tolerability. Results. The post-titration cohort comprised 24,876 confirmed slow responders. Within-class dose escalation produced a small consistent benefit at 24 weeks: AIPW +0.64 pp (95% CI +0.16 to +1.12), with five non-AIPW estimators ranging +0.47 to +0.76 pp. The continue arm itself lost an additional 8.07 pp over the same window (96% continued to lose), so escalation is a marginal addition to a substantial natural slope, not a rescue. Equipotent class switching from semaglutide to tirzepatide was inconclusive: linear and matching estimators ranged +1.26 to +1.80 pp, but AIPW was -0.33 pp (95% CI -1.27 to +0.60) with only 90 treated patients and limited propensity-score overlap. Class switch with simultaneous potency increase (sema to tirz) gave AIPW +0.65 pp at 12 weeks (95% CI +0.37 to +0.92, n = 80). A negative-control specification yielded ATE -0.13 pp, indicating the pipeline did not generate spurious signal. A held-out-fold prognostic-threshold sensitivity gave a null effect (+0.04 pp), correcting an earlier circular +1.06 pp estimate. Conclusions. Among confirmed post-titration slow responders, within-class dose escalation adds approximately 0.6 percentage points at 24 weeks on top of an 8 percentage point natural slope, consistently across six estimators including doubly-robust inference. This headline effect is small and not robust to modest unmeasured confounding (E-value 1.27) or to MNAR-style outcome attrition (tipping point delta approximately 1.2 pp), so it should be read as hypothesis-generating rather than practice-changing. Class-switching evidence is inconclusive; linear-estimator results suggesting benefit did not survive doubly-robust estimation in small treated samples with limited propensity overlap. The dose-ladder framework, with phase-specific evidence grading, is hypothesis-generating and insufficient on its own to change practice.

Background: Urinary tract infection (UTI) is the most common secondary condition among people with spinal cord injury/disease (SCI/D). Intravesical Lacticaseibacillus rhamnosus GG (LGG) is an antibiotic-sparing approach to managing urinary symptoms. Objective: Determine the optimal number of doses of intravesical LGG for urinary symptom reduction. Design: Prospective, randomized, two-arm dosing trial. Setting: National recruitment with a local subsample providing urine samples in Washington, DC, USA. Participants: Adults with SCI/D and neurogenic lower urinary tract dysfunction (NLUTD) who use intermittent catheterization (IC); 177 enrolled and randomized (intention-to-treat), with 76 compliant instillers (39 low-dose, 37 high-dose) in the per-protocol analytic sample. Interventions: Two (2 doses/24 hours) or four (4 doses/36 hours) intravesical LGG regimens, self-initiated in response to cloudier or malodorous urine per the Self-Management Protocol using Probiotics (SMP-Pro). Main Outcome Measures: Primary: proportion achieving [&ge;]20% reduction on the Urinary Symptom Questionnaire for Neurogenic Bladder-Intermittent Catheter version (USQNB-IC). Secondary: urinary biomarkers (leukocyte esterase, nitrite, white blood cells, urinary neutrophil gelatinase-associated lipocalin [uNGAL]) and standard urine culture (SUC) in a local subsample. Results: By Day 2, 57.9% (63.8% low-dose; 51.2% high-dose) achieved [&ge;]20% total symptom reduction; high-dose success rose to 70.0% by Day 4. Thirty percent of high-dose participants did not respond at either time point and could not be distinguished from responders by demographics or urine biomarkers. Urinary biomarkers and SUC were unchanged pre- to post-instillation. No serious adverse events were adjudicated as attributable to intravesical LGG by an independent Data Safety Monitoring Board (DSMB). Conclusions: A two-dose course of intravesical LGG yields clinically meaningful symptom improvement in the majority of people with SCI/D and NLUTD who use IC; four doses benefits a meaningful subgroup of two-day non-responders, while a small cohort remains nonresponsive. These results provide preliminary dosing guidance and support progression to a definitive trial.

Objective: There is a bidirectional relationship between emotional distress and pain. However, this relationship is understudied in people with HIV in low-resource settings. This study sought to describe the temporal relationship between emotional distress and pain in people with HIV. Design: Longitudinal observational study. Methods: Participants with virally suppressed HIV, reporting either no pain or persistent pain at baseline, provided weekly remote ratings of distress, worst pain, and average pain using 0-10 visual analogue scales. Within-individual fluctuations in distress and pain were visualised over time. Group-level correlations were determined using Spearman's correlation tests. Cumulative link mixed models assessed whether distress and pain each predicted the other in the following week. Results: 72 participants provided responses over 49 weeks. The participants had a median (IQR) age of 43 (37-51) years, 63% (n=45) were unemployed and most were females (n=51;71%). Distress and pain fluctuated concurrently within individuals: distress was positively correlated with worst pain ({rho}=0.66, 95% CI= 0.60-0.72, p<0.001) and average pain ({rho}=0.70, 95% CI=0.64-0.75, p<0.001) intensity within the same week. Worst pain (OR=1.42, 95% CI=1.17-1.71, p<0.001) and average pain (OR=1.43, 95% CI=1.20-1.71, p<0.001) intensity both predicted distress in the next week. Distress predicted worst pain intensity (OR=1.25, 95% CI=1.07-1.46, p=0.023) but not average pain intensity (OR=1.19, 95% CI=1.01-1.40, p=0.152) in the next week. Conclusions: The temporal relationship between distress and worst pain intensity was bidirectional, whereas distress did not temporally predict average pain intensity. Both pain and emotional distress should receive attention from HIV research and clinical care in low-resource settings.

Background: Actigraphy-derived rest-activity rhythm (RAR) features are widely used to characterize clinical states in bipolar disorder (BD). Both mean levels and temporal variability of these features have been associated with mood episodes; however, variability measures are often statistically coupled with the mean, particularly in skewed distributions. This raises a question as to whether variability reflects a separate characteristic of the data or whether the observed association arises from statistical properties of the data. Objective: In this study, we aim to determine whether temporal variability of actigraphy-derived RAR features provides standalone information on mood episodes in BD beyond mean activity levels after accounting for mean-variance dependence. Methods: We analyzed actigraphy data from a subset of 72 participants with BD drawn from a larger longitudinal study, extracting 22 daily RAR features aggregated weekly as sample mean (MEAN) and within-week temporal variability computed as sample standard deviation (VAR). Variance-stabilizing transformations (Box-Cox or Yeo-Johnson) were applied to the entire study cohort to reduce mean-variance dependence. Associations with mood episodes and remission (mania: n=34; depression: n=58 annotated participants) were evaluated using generalized linear mixed-effects models with a logistic link function, including univariate (MEAN or VAR) and multivariate (MEAN+VAR) specifications, assessed by likelihood-based metrics and the area under the receiver operating characteristic curve (AUC). Results: Transformations reduced mean-absolute correlations from 0.43 to below 0.06. Temporal variability remained significantly associated with clinical state for 11/22 RAR features in mania and 16/22 features in depression, with all significant associations remaining after false discovery rate correction (p<0.05). Joint models showed modest incremental gains (AUC 3%-4% overall; up to 12% in mania, 7% in depression), with absolute performance remaining limited (AUC 0.50-0.66). In both mania and depression, nearly all significant variability-based regressors contributed incremental information beyond mean-based models. Only sleep duration and activity changes around wake time (+-1 hour), did not improve discrimination between mania and remission. Conclusions: Temporal variability in RAR features can be considered a standalone state marker of mood episodes not captured by mean activity. We found it to be more consistently associated with depression than mania. Its incremental discriminative contribution is modest, suggesting greater utility within multivariate or multimodal frameworks.

Background: Children born prematurely (before 37 weeks) or admitted to the neonatal unit (NNU) are at increased risk of adverse long-term physical health outcomes. It is also recognised that there is an association with later academic performance and special educational needs, however it is not clear whether these broad risk factors could be used as stand-alone heuristics to identify children who may benefit from additional support in educational settings. We aimed to examine the associations between neonatal unit (NNU) admission and educational attainment in mid-childhood. Methods and Findings: Pregnancy data from a prospective birth cohort (Pregnancy Outcome Prediction Study, Cambridge, United Kingdom, 2008-2012) were linked to national educational outcomes (Department for Education, United Kingdom). Multivariable regression models adjusted for maternal, child, and socioeconomic factors were used to evaluate associations between (i) all NNU admissions, (ii) at term NNU admissions >48 hours, (iii) preterm birth without ongoing physical health needs, and educational outcomes at ages 5-11 years. Children who required any NNU care were more likely not to meet expected educational standards across multiple ages and domains in early and mid-childhood: age 5 early year foundation (aOR 1.64, 95% CI 1.19-2.27, p=0.003), phonics at age 6 (aOR 2.43, 95% CI 1.72-3.57, p<0.001), and at age 7 (here assessments were divided into multiple domains): reading (aOR 1.67, 95% CI 1.18-2.38, p=0.004), writing (aOR 1.72, 95% CI 1.25-2.38, p<0.001), mathematics (aOR 1.56, 95% CI 1.09-2.22, p=0.020), and science (aOR 1.85, 95% CI 1.22-2.78, p=0.003). Similar patterns were observed among both at term-born infants who stayed >48hrs in NNU (phonics assessment at age 6 aOR 2.26, 95% CI 1.51-3.36, p<0.001) and in children born preterm without long-term physical health sequelae (phonics assessment at age 6 aOR 3.07, 95% CI 1.96-4.81, p<0.001). These associations were robust to adjustment for demographic, perinatal, and socio-economic factors. By age 11, differences in academic attainment were attenuated and no longer clearly distinguishable across all exposure groups. However, there was an increased likelihood of special educational needs (SEN) at age 11 associated with any NNU admission (aOR 1.78, 95% CI 1.15-2.73, p=0.009), at term NNU admission for >48hrs (aOR 1.88, 95% CI 1.19-3.00, p=0.007), and children born preterm without long-term physical health sequelae (aOR 1.50, 95% CI 1.00-2.25, p=0.049). Predictive performance of any NNU admission for SEN at age 11 was moderate (AUC 0.70, 95% CI: 1.14-2.65, p=0.010), with balanced sensitivity and specificity and high negative predictive value. Conclusions: NNU admission, for both term and preterm infants, is associated with poorer educational outcomes and an increased likelihood of special educational needs in mid-childhood.

Background The 2022-2023 global mpox outbreak highlighted the critical need for robust genomic surveillance capabilities to track mpox virus (MPXV) evolution and transmission dynamics. Methods Building upon our established SARS-CoV-2 sequencing infrastructure, we implemented a Molecular Loop probe-based long-read sequencing approach using Pacific Biosciences Sequel II technology for comprehensive MPXV genomic surveillance across the United States (US). From August 2024 to June 2025, we generated 326 high-quality whole genome sequences from residual mpox-positive clinical specimens collected by Labcorp across all 10 US Department of Health and Human Services regions. Results Our analysis identified two samples containing clade Ib MPXV in January and June 2025 and captured shifting trends in clade IIb diversity, with 13 distinct lineages observed. We also identified multiple instances of large (~1.6-17.6kb) deletions proximal to the inverted terminal repeats in clade IIb genomes. APOBEC3 mutation analysis indicated substantial evidence of human-to-human transmission among both clades. Further, we observed significantly higher APOBEC3-associated SNPs per kilobase (P<0.001) in clade IIb genomic variable regions relative to their central conserved region. Our assay exhibited strong reproducibility across biological replicates from individual patients and accuracy was confirmed via parallel sequencing of select specimens by US Centers for Disease Control and Prevention (CDC) using metagenomic sequencing. We also demonstrated via custom simulation that our assay discriminates all known MPXV clades and lineages, including those we have not observed in the US. Conclusions Our integrated nationwide surveillance system facilitates real-time genomic tracking of outbreak evolution, with demonstrated capacity across SARS-CoV-2 and MPXV, positioning this platform for rapid deployment during future pathogen emergence.

Background: Monoclonal antibodies have emerged as a prophylactic strategy to prevent symptomatic SARS-CoV-2 infection in immunocompromised individuals. However, the evolutionary and clinical implications of breakthrough infections under this regime remain unclear. Methods: A male in their 80s with a haematological/oncological diagnosis received a 2000 mg intravenous infusion of sotrovimab in March 2023 and was diagnosed with COVID-19 by RT-qPCR from a nasopharyngeal swab in August 2023. Weekly samples (n=24) were collected through February 2024 (171 days). All samples underwent whole-genome sequencing, with select mutations subjected to functional assessment. Findings: Sequencing identified the GE.1 lineage at all timepoints. An intra-host recombination event in ORF1ab (positions 8942-12458) was detected prior to 23 weeks post-detection, followed by a 14-fold increase in viral load (7.42e+06 to 1.00e+08 RNA copies/mL) and a marked shift in the viral population. E340D, a sotrovimab resistance mutation, was detected at low abundance (46%) within the first week post-infection, fluctuated over time, and was nearly fixed by week 15 (107 days) post-detection. We assessed five spike mutations - V36M, S98F, and V213G in the N-terminal domain, Y505P in the receptor-binding domain, and P681Q near the S1/S2 cleavage site - and additionally evaluated the impact of E340D. V36M conferred the highest infectivity across all cell lines, with the most significant effect in low-TMPRSS2 cells. While all mutations showed enhanced infectivity with the addition of E340D, the effect was most pronounced in mutations with lower baseline infectivity. The addition of E340D significantly decreased relative neutralizing titres for V36M, S98F, and V213G, enabling escape from neutralizing antibodies in XBB-responsive individuals, illustrating an enhanced phenotypic advantage. Patient neutralizing activity was absent pre-sotrovimab, and sotrovimab-induced neutralization was further compromised by selection of E340D. Interpretation: Sotrovimab pre-exposure prophylaxis in an immunocompromised patient did not prevent SARS-CoV-2 infection, and selected for resistant mutation E340D, with unexpected fitness consequences across non-receptor binding domain spike regions.

Background: Opioid use disorder (OUD) is heritable, yet most genome-wide association studies (GWAS) have focused on European populations, leaving the genetic architecture of OUD in non-European populations underexplored. Methods: We conducted GWAS of OUD across three ancestries using electronic health records and genomic data from 52,357 All of Us Research Program participants (8,912 cases; 43,445 matched opioid-exposed controls; 48.5% female). Participants were stratified into European (EUR), African (AFR), and Admixed American (AMR) ancestry groups for logistic regression GWAS, with independent replication in the Million Veteran Program. We then applied the deep-learning model AlphaGenome to predict the tissue-specific transcriptomic and splicing consequences of top risk variants across 13 reward-pathway brain regions. Results: We identified and replicated a novel DDX6 risk locus, alongside established OPRM1 and FURIN signals. AlphaGenome predicted the DDX6 regulatory allele downregulates the stress-resistance gene FOXR1 in the nucleus accumbens, while the protective OPRM1 variant (rs1799971) upregulates OPRM1 expression across reward networks. Other signals of interest included IL6R and SHISA9 (EUR); GHR (AFR); and ASTN2 (AMR). Conclusions: This study identifies DDX6 as a novel OUD risk locus, replicates associations with OPRM1 and FURIN, and highlights biologically plausible ancestry-specific signals in AFR and AMR populations. We also replicated top variants in an independent population. Finally, integrating GWAS with deep-learning annotations provides specific, localized biological hypotheses to guide future experimental validation and targeted therapeutics.

Stochastic epidemic models are a cornerstone of infectious disease epidemiology and are often used to study intervention scenarios. However, large run-to-run variability can make intervention effects difficult to estimate precisely. We revisit the epidemic Sellke construction, which assigns each individual an infection threshold for the cumulative infection hazard such that, conditional on the thresholds, the epidemic trajectory becomes deterministic. This enables coupling of simulations with and without an intervention, yielding low-variance effect estimates even when outcomes such as final size or peak incidence vary widely between runs. We develop an exact, event-driven implementation that maintains infection and recovery events in priority queues. Cumulative infection-hazard updates require O(log N) time per event, yielding overall complexity O(Elog N) for E events in a population of size N. The implementation achieves computational performance comparable to the classical Gillespie algorithm while naturally accommodating non-Markovian infectious periods and complex infectiousness profiles. We illustrate the approach using distance-dependent spread of avian influenza between poultry farms in the Netherlands and a multilayer population with households, schools, and workplaces. In both examples, coupling enables efficient within-run comparisons of intervention scenarios across stochastic realisations.

Estimating dengue force of infection (FOI) is essential for understanding transmission dynamics and targeting intervention programmes, yet surveillance data in endemic settings required for estimations are often incomplete, with varying formats. We developed a Bayesian hierarchical catalytic model that jointly fits age-stratified case data, aggregate case data, and seroprevalence surveys within a single framework, incorporating external covariates to improve parameter identifiability. Synthetic validation showed that covariates alone recovered accurate FOI point estimates even when most districts contributed only aggregate data, but did so with poorly calibrated uncertainty; anchoring the model with a single seroprevalence survey was necessary to bring credible interval coverage close to nominal. Applied to 128 districts across Java and Bali, Indonesia (2016-2024), the model revealed substantial spatial heterogeneity in FOI and reporting rates. Many districts in Java exceeded the WHO-suggested seroprevalence threshold for vaccine introduction, yet were classified as low-priority when using reported incidence as prioritisation criterion, particularly in areas with weak surveillance. Model-based seroprevalence estimation, integrating multiple data sources, offers a more consistent basis for identifying high-priority districts for vaccine introduction, and is less susceptible to surveillance bias than reported incidence.

Background Post-infection conditions (PICs), such as Long Covid, are associated with heterogeneous, fluctuating symptoms that profoundly affect daily functioning. Despite moderate-certainty evidence from the NIHR-funded LISTEN trial (COV-LT2-0009) that personalised self management support improves outcomes and may reduce societal and economic impacts of Long Covid, many people living with PICs still receive condition-specific services, generic advice, or stand-alone digital tools that do not address their complex needs. Aim To map care approaches in general practice and synthesise UK evidence for PIC management. Design and setting Scoping review and online survey. Method A two-phase study was conducted: (1) a scoping review of UK evidence on PIC management in general practice; and (2) a supplementary online survey of practitioners working in UK general practice to provide contextual insights. Results The scoping review identified 32 studies focused on Long Covid. One study included a comparator group (ME/CFS). Study populations were predominantly white ethnicity and female. Evidence for non-Covid PICs in UK general practice was largely absent. The supplementary survey (n=46) provided preliminary practice-level insights. Healthcare practitioners reported varied PIC presentations, diagnostic uncertainty, limited referral pathways, inequitable access, and low confidence in managing PICs. Conclusion Evidence informing PIC management in UK general practice remains predominantly Long Covid-focused and may not reflect the range of PICs encountered in practice. While survey findings are preliminary and require confirmation in larger samples, they highlight uncertainty around PIC management. Further research is needed to evaluate whether existing Long Covid pathways should be expanded or complemented by broader PIC models. Keywords general practice; Long Covid; self-management; post-viral syndromes

Introduction: The Diabetes Prevention Program Outcomes Study (DPPOS) is an established cohort of aging persons with pre-diabetes and type 2 diabetes with 25 years of median follow-up. In 2022 DPPOS added Alzheimer's disease (AD), and AD related dementias (ADRD) phenotyping using the National Alzheimer's Coordinating Center (NACC) Uniform Data Set (UDSv3), which included a standardized neurological examination across 25 clinical sites, administered by clinical staff and interpreted centrally by clinicians. Methods: A DPPOS video-based asynchronous neurological examination (DPPOS-VANE) was developed iteratively through consensus from research clinicians and staff feedback to harmonize with UDSv3 to identify common neurological diagnoses aside from dementia including diabetic cranial neuropathies, stroke and parkinsonism. DPPOS-VANE was designed to be conducted without direct participant contact by the examiner, reproducible, and independent of clinical skills of PCs. An iPad camera recorded the video exam, comprised of assessments of extraocular and facial movements, visual fields, speech, gross motor strength, pronator drift, praxis and parkinsonism. A 10-minute training video demonstrated the examination step-by-step with scripts and instructions in English and Spanish. Site-specific performance review, feedback, and staff certification preceded central reading of video recordings by physicians. After two years of implementation, 1286 DPPOS-VANEs led to 1284 examination reviews. Of these, 1204 (93%) were completed by having the examiner follow the standard script. Overall, 1237 examinations (96%) were delivered as planned, 41 (3%) had minor errors but were still usable, and 6 (0.4%) had major deviations in exam technique; two additional recorded evaluations were not usable as recorded videos were inaccessible due to technical errors. Each examination was completed within 10-15 minutes. Each site on average completed 51.4 examinations (range 14-92). Discussion: Engaging 55 research staff across 25 sites and 3 physician-reviewers, this study is the first to demonstrate feasibility of a VANE as an efficient neurological examination model enabled by commonly used devices. Such a multisite standardized VANE represents a novel paradigm for large epidemiological studies.

Controlled human infection models (CHIMs) are increasingly used in infectious disease research to study pathogen dynamics and evaluate interventions under controlled conditions. However, these studies are resource-intensive and involve ethical and safety constraints, making efficient study design critical. Dose-finding is a key early component in CHIMs, where the aim is to identify a challenge dose that achieves a target infection probability. Traditional rule-based designs are commonly used but can be inefficient, motivating the use of model-based adaptive approaches such as the Bayesian Continual Reassessment Method (CRM). Although CRM has been extensively studied and widely adopted in Phase I oncology trials for identifying the maximum tolerated dose of therapeutics, its application in CHIM settings remains limited, particularly when the endpoint of interest is infection. This tutorial provides step-by-step guidance for implementing a Bayesian CRM in dose-finding CHIMs, using an oropharyngeal Neisseria gonorrhoeae challenge as a motivating case study. The framework outlines key design components, including dose-grid specification, dose-response model, prior elicitation, Bayesian updating, decision rules, and stopping criteria, with particular emphasis on a clinically interpretable parameterisation. Trial operating characteristics are evaluated through simulation studies under multiple dose-response scenarios and prior-predictive analyses, and compared with a commonly used '3+3' type rule-based design. This work highlights the advantages of Bayesian model-based designs for dose-finding in CHIMs over classic rule-based designs and provides a structured, reproducible framework for implementing CRM, supporting their application in future CHIM studies.

Lithium has been the gold standard for the treatment and prevention of relapse in bipolar disorder for over 60 years. Guidance from the National Institute for Health and Clinical Excellence states explicitly to 'offer lithium as a first-line, long-term pharmacological treatment for bipolar disorder'. Yet, in the last two decades its use has been in decline with clinicians favouring anticonvulsants or antipsychotics when treating this condition. In this study, we have used three openly available datasets containing prescribing data from primary and secondary care to explore trends in the use of lithium in England, showing both regional and temporal variance between 2015-2024. We have shown that lithium use declined in primary care by 20.9% in the last ten years (2015-2024) and 10.9% overall in the last five years (2019 to 2025). We have also shown how there is some regional variation in the source of lithium for patients, although the vast majority is prescribed in primary care. Further research into clinical behaviour is needed to understand what is driving the decrease in lithium usage, and what barriers and enablers may influence its use across the country.